5 



Oh! 



The global colour model of QCD and its relationship to the 
NJL model, chiral perturbation theory and other models * 

Reginald T. Cahill and Susan M. Gunner 
Department of Physics, Flinders University ^ 

GPO Box 2100, Adelaide 5001, Australia 
and Institute of Theoretical Physics, Adelaide 



a^ ■ hep-ph/9601319 

O 
(N 

(N 
> 
0^ 



Abstract 



The Global Colour Model (GCM) of QCD is a very successful model. 
Not only is it formally derivable from QCD but under various conditions 
i it reduces to the NJL model and also to Chiral Perturbation Theory, and 

^ \ to other models. Results presented include the effective gluon propaga- 

tor, the difference between constituent and exact quark propagators, 
various meson and nucleon observables, a new form for the mass for- 
mula for the Nambu-Goldstone mesons of QCD, and the change in the 
MIT bag constant in nuclei. 



1. The Global Colour Model (GCM) 

The GCM can be formally derived from QCD Q. The GCM is in turn easily related to 
^ I a number of the more phenomenological models of QCD as indicated in fig.l. 
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QCD and models for QCD amount to the calculation of correlation functions. 



G{.-,x,...) = J T>qT>q'DA....q{x) exp{-Sqcd[A,q,q]) (1) 

At low energies we expect that only hadronic correlations are directly observable, and we 
expect (1) to be transformable into an hadronic form, equivalent to the effective action 
description of low energy nuclear physics. 



j VqDqVAexp{-Sqcd[A, q, q] + rjq + qr]) 

J VTTVNVN...expi-Shad[TT,...,N,N,..] + J^[ri,ri]n + ..) (2) 
If in (1) we formally do the gluon integrations, after fixing a gauge (ghosts not shown) 

J VqVqVAexp{-Sqcd[A,q, q]+r]q + qr]) = j 'DqVqexp{-S[q, q]+r]q + qr]), (3) 
leaving an action for quarks of the form 

S[q, q]=l qil.d + M)q + \j f,ix)fAy)D,Ax - v) + f.jtfpDfj^p + - (4) 
where = q{x)^jf^q{x), and D^jj{x) is the exact pure gluon propagator 

D^,{x) = j VAAlAlexp{-Sq,d[AM) (5) 

A variety of techniques for computing D^jy{x) exist, such as Dyson-Schwinger equations 
(DSE) lattice simulations Q and, as discussed in other papers at the workshop, 
monopole and instanton modellings. Dropping the D'^p, .. defines the GCM. This is 
equivalent to using a quark- gluon field theory with the action. 

Sg,m[q,q,A';,] = J i^q{j.d + M+zAl^^,)q+^A';^D-^{zd)A:^ (6) 

where D'^^{p) is the matrix inverse of D^y{p). This action has a global colour symmetry, 
hence the description as the GCM. Hadronisation of the GCM involves a sequence of 
functional integral calculus changes of variables involving, in part, the transformation to 
bilocal meson and diquark fields, and then to the usual local meson and baryon fields. 



VqVqVAexp{-Sqcd[A,q,q]) (7) 

« j VqVqexp{-J q{l.d + M)q + j jl{x)jl{y)D^,{x - y)) (8) 

= J DBDVDV*exp{-S[B,V,V*]) (bilocal fields) (9) 

= JvTrVNVN...exp{-Shad[Tr,-,N,N,..]) (local fields) (10) 
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The derived hadronic action, to low order in fields and derivatives, has the form 



Shad[Tr,...,N,N,..] = 

J (fxtr{N{-f.d + mN + Atun - mNyj2i-i^T:°-T'' + ..)N} + 



+ / (fx 
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ifulfpflTGuipTY€fj_i,o-T^fJ.di>Po-'dT-TT~\~ 

(11) 



+ ^e^yartr{TT.Fd^TT.Fdy7r.Fd^7r.Fdr7T.F) + .... 

This induced effective action is the action of Quantum Hadro-Dynamics (QHD). Being a 
derivative expansion it should not be used in hadronic loop calculations. For that purpose 
the non-local form of the hadronic effective action is necessary. The above derivation 
leads naturally to the dynamical breaking of chiral symmetry (sect. 2) and so to the Chiral 
Perturbation Theory phenomenology (ChPT) (sect. 4). The bare mesons are described by 
Bethe-Salpeter equations (BSE), and the baryons by covariant Faddeev equations (FE) 
Q. The GCM thus relates hadronic properties directly to the gluon D^y. It is not known 
why the GCM truncation is so effective. 

2. Dynamical Chiral Symmetry Breaking 

The hadronic effective action in (10) arises from expanding the bifocal action in (9) 
about its minimum via the Euler-Lagrange equations (ELE): 6S/6B = give the rainbow 
constituent quark DSE in (12) and (13), while 5S/5T> = has solution T> = 0. Fluctu- 
ations are described by the curvatures: 5'^S/6B6B gives the ladder BSE for mesons, and 
6'^S/6V6V* gives the ladder BSE for diquarks. 

, 16 f d^q , B(q)+m 

'■^fe'> - = 1 1 1?^-""'" - ' \^Ai,r fit) + mr 

G{q) = iiA{q)q.-f + B{q) + m)"! = -iq.-faM) + a^iq) (14) 

Note that the quark propagator in (14) is a constituent quark propagator (and appears 
in the exponent in (9)), and not the exact quark propagator, as would arise from (1), 
and which would satisfy the Slavnov- Taylor identities (STI). In the GCM and related 
phenomenologies the quarks in the bare hadrons in (11) are further dressed, for example 
by mesons, in the functional integrations of (10). Imposing the STI in (12) and (13) 
could thus result in double counting problems. The exact quark propagator is relevant to 
the question of absolute quark confinement. The constituent quark propagator has the 
running mass M{q) = B{q) / A{q). The constituent quark mass arises as the most likely 
value of M{q) in, say, the BSE. It will depend on the bare current mass and to some extent 
even the hadron state. 
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3. The Nambu - Jona-Lasinio Model (NJL) 

The NJL model is a special case of the GCM. Formally the NJL model |^ is the contact 
interaction limit of the GCM: D^y{x — y) ^ g'^5^y5{x — y) or D{p) g^. But in the 
constituent quark DSE the contact limit is undefined because it leads to divergences in (12) 
and (13). A cutoff A is then always introduced in NJL computations, which is equivalent 
to using the gluon propagator D{p) g^Oiq^ — A^). Hence the NJL model is the GCM 
but with a box-shaped D{p), rather than a 'running' D{p). The GCM bosonisation in (8) 
arises from a generalised Fierz identity 

= -9(^)^^(x)^^g(x)^^^g(x) (15) 

where {M^} = W^K"" O 1 F^} and {M^} = { Vf^"' O i?^} involve generators of 
the Dirac, colour and flavour symmetries. The use of this identity gives the generalized 
NJL model. 

4. Chiral Perturbation Theory (ChPT) 

When the quark current masses M. ^ S[q,q,A'^] has an additional global Ul{Np) ® 
Ur{Nf) chiral symmetry: writing 97^9 = qB^^qR+qLl^l1L where qR^i = PrM and qj^ ^ = 
qPl,r we see that these two parts are separately invariant under qn — > URqR,q^ — > qRUj^ 
and qL — > UiqLjQL ~^ ^L^l- consequences may be explicitly traced through the 
GCM hadronisation. First the ELE 6S/6B=0 have degenerate solutions. In terms of the 
constituent quark propagator we find 

G{q; V) = [iA{q)q.j + VB{q)r' = C^G{q; 1)C^ (16) 
where 

C = VV, y = exp(V275^"i^") (17) 

and consequently the fluctuations 6'^S/5BSB have massless Nambu-Goldstone (NG) BSE 
states. 

In the hadronisation (10) new variables are forced upon us to describe the degenerate 
minima (vacuum manifold): 

U{x) = ea;p(iV2^"(x)F") (18) 
V{x) = PlU{x)^ + PrU{x) = exp(V275vr"(x)F'^) (19) 
The NG part of the hadronisation gives 



+ 



+K2tr{[df,ud^u^f) + KMd^ud,u^d^ud,u^) + ....) (20) 



This is the ChPT effective action Q, but with the added insight that all coefficients 
are given by explicit and convergent integrals in terms of A and B, which are in turn 
determined by -D^jy. The higher order terms contribute to vrvr scattering. The dependence 
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of the ChPT coefficients upon D^j^u have been studied in ||^, The hadronisation 

procedure gives a full account of NG-meson - nucleon coupling. 

5. The MIT, Cloudy Bag and Soliton Models 

While the GCM hadronisation in (11) and (20) is the main result, at an intermediate 
stage one obtains [|l|] extended meson Quark- Meson Coupling type models (QMC) 
Applying mean field techniques to this GCM quark-meson coupling effective action leads 
to soliton type models, which have been studied in detail in [^] and the significance of the 
extended mesons demonstrated. From the soliton models a futher ansatz for the soliton 
form leads to the MIT and Cloudy Bag Model (CBM). In particular the GCM expression 
for the MIT Bag constant may be extracted: 

^ 127r2 , A'^(s)s + B^(s), B^(s) , , , 

S = / sds[ln( — ^ 'J^ — ^) - , ^ J (21) 

(27r)4yo A2(s)s ^ A^{s)s + B^{s)^ ^ ' 

which is based on the energy density for complete restoration of chiral symmetry inside a 
cavity. This bag constant is for core states as no meson cloud effect is included. 

With a mean field description of the pion sector via (t{x)^ which describes the isoscalar 
part of g{x)V{x\ where ct(x) is a 'radial' field multiplying the NG boson field V{x) 

, 127r2 /-^ , A^{s)s^a^B'^[s)^ a^B^{s) ^ , , 

^(^) = (2^yo ^^^[^^^T(^)-lW^^ ^''^ 

which reduces to (21) when a = 1, being the non-perturbative field external to an isolated 
nucleon core, and a < 1 describing a partial restoration of chiral symmetry outside of 
the core. Using the gluon propagator discussed in sect. 6 we obtain the plot of B{a)/B{0) 
shown in fig. 2. Dressing of the nucleon core by mesons is partly described by a reduction 
in a in the surface region, causing a reduction in the nucleon mass. 
However in nuclei a mean meson field description [14| means that a is even further 



reduced outside of the nucleons, and the effective bag constant is further reduced. The a 
field, which can model in part correlated vrvr exchanges, along with the u; meson field, are 



believed to be important to a mean field modelling of nuclei. In [15| it has been argued 
that the reduction of the effective bag constant for nucleons inside nuclei is essential to 
the recovery of features of relativistic nuclear phenomenology. The GCM thus allows B{a) 
and details of relativistic nuclear phenomenology to be directly related to the constituent 
quark propagator, and in turn to the gluon propagator. 

6. Gluon Propagator Separable Expansions 

In application the GCM model amounts to solving a sequence of non-linear and linear 
integral equations: 
D^u Gconst.quark ^ meson BSE & diquark BSE — > 
covariant Faddeev eqns for baryons hadronic observables 

In GCM computations one uses what amounts to a mixed metric: a Euclidean metric 
for the internal quark-gluon fluctuations, and with the hadron momenta in the time-like 
region of the Minkowski metric (usually but not necessarily in the centre-of-mass (cm) 
frame). The constituent quark propagator equations (12) and (13) are solved in the 
Euclidean region {s = q"^ > 0), and then A{q'^) and B{q^) are 0(4) invariant. Significantly 



it was discovered [16| that when low mass BSE states are solved with this mixed metric 
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approach the (on-mass-shen) bound state form factors also show a remarkable degree of 
0(4) invariance wrt the relative quark momentum, even though this was not assumed in 
the numerical solution of the BSE. If we note that an 0(4) hyperspherical expansion 

D{p -q)= Doip^ q^) + q.pDi{p^, q^) + ... (23) 

of the gluon propagator, where 

Dq{p^, q^) = - r dp sin'^p D{p^ + q^ - 2pqcosP) (24) 
vr Jo 

together with the numerical technique of using a multi-rank separable expansion 
Do{p\q')= E Up')Uq'), (25) 

i=l,n 

then this separable expansion automatically generates 0(4) invariant BSE solutions when 
(23) is used, at lowest order, in a BSE. Hence the significant realisation that using a multi- 
rank separable expansion is ideally suited for GCM computations. This technique then 
renders the DSE and BSE equations to essentially algebraic form, with only the baryon 
integral equation computations requiring extensive numerical solution. As shown in [^] 
this very useful property of 0(4) BSE solution invariance appears to be a consequence of 
having no nearby singularities in the constituent quark propagator, and so appears to be 
related to the confinement property of the quarks. 

It is important to note that the separable technique is introduced only as a numerical 
technique and only after all integral equations have been formally derived using the co- 
variance of the GCM. If the separable expansion of the gluon propagator were introduced 
in the defining action of the GCM then the explicit breaking of covariance would block 
any of the usual momentum-space DSE, BSE,., equations from being derivable. 

Implementing the separable expansion we first solve the DSE (12) and (13). Then 
sums are obtained 

B{s) =^biTi{s), ... asis) = ^ as{s)i, cr^is) = ^ ay{s)i (26) 

i=l,n i=l,k 

The gluon propagator is implicitly parametrised in |^ by assuming entire functions (as a 
means of implementing quark confinement) for ag and a^, say 

/ , N 2s- P^(l- exp(-2s/p'^)) 
(7s{s)i = Ciexp{-diS), o■^,(s) = (27) 

Then (12) gives 

.2 16_2 B(,,) 



in which B{s) = B{s)i + -6(5)2 + ••, and B{s)i = as{s)i/{sa.a{sy + ^^(s)^). 

Having solved the constituent quark non-linear equations (12) and (13) one can then pro- 
ceed to solve the meson and diquark BSE and finally the constituent nucleon-core Faddeev 
equations |l|]. The implicit gluon propagator parameters in (27) are then determined by 
a best fit to some hadron data: fT^,m-,^,ma^. Subtleties associated with the incorporation 
of the quark current masses and the reasons for choosing to fit these particular data are 
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discussed in |8|. The gluon propagator parameter values are also given in and various 
hadron observables are presented, and discussed here in sect. 7. 

Having determined the best fit the translation invariant form of the gluon propagator 
may be determined from (24) which, with = 0, gives 



z^(p2) = Do(p^o) = 5]^,(p2)^,(o) = 

i 

However the form of the extracted gluon propagator will clearly depend on the forms 
assumed in (27), and futhermore, while a unique best fit was obtained the fit is somewhat 
flat wrt variation of some of the parameters. To help more accurately resolve the form of 
the effective gluon propagator that arises in the GCM, current work |17] that constitutes a 



Hybrid Lattice-GCM calculation. In this a multi-rank separable expansion is constructed 
and its parameters determined by fits to both some hadron data and as well to the form for 
the gluon propagator that arises from lattice computations except for the infrared 

(IR) region of the gluon propagator. This amounts to probing the deep IR properties of 
the gluon propagator by means of low energy hadronic data. Results will be reported 
elsewhere. 

7. Meson and Nucleon Observables 

Table 1 shows various hadronic observables computed using the separable expansion 
technique of sect. 6. Of particular relevance here are the tttt scattering lengths which 
arise from the ChPT effective action (20) with the parameters given by the explicit GCM 
forms 0. The constituent quark masses arise from the value of the constituent quark 
running mass M(s) at the value corresponding to the dominant s value in the BSE. This 
can only be determined after the BSE solution is known. Various diquark masses are 
shown. Diquarks are extended quark-quark correlations in baryons, and arise naturally 
from the GCM hadronisation They are particularly effective in reducing the numerical 
complexity when solving the baryon Faddeev integral equations. Note that a nucleon-core 
mass of 1390MeV is consistent with the first Faddeev computation of the nucleon-core 
mass in [^]. This is also consistent with the idea that dressing of the nucleon-core by 
NG bosons causes a decrease of some SOOMeV. Incorporation of the spin 1+ diquark state 
into the nucleon-core computation is also expected to decrease the nucleon mass. Other 



computations of the nucleon-core state [19, 21, 22 1 in the context of the NJL limit of 
the GCM have always adjusted the NJL parameters so that the nucleon-core mass actually 
fitted the experimental nucleon mass, leaving out the significant NG dressing effect. 

Also shown is the GCM predicted value for the MIT bag constant, using expression (21). 
This is somewhat larger then the usual MIT value and gives a MIT nucleon-core mass of 
1500MeV without cm corrections. With cm corrections we might expect this to reduce 
to near the covariant Faddeev value of 1390MeV, i.e the GCM MIT bag constant value 
appears to be consistent with the expected nucleon-core mass. The further reduction in the 
GCM B value when chiral symmetry is partially restored due to meson dressing of isolated 
nucleons, and a further enhanced reduction for nucleons inside nuclei, was discussed in 
sect. 5. 
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8. A New NG Mass Formula 

The usual NG mass formula is 

where 

p =< qq >= NMG{x = 0)) = U J ^^s(g') (31) 

This expression for p is divergent in QCD. The values of m and <qq> are then usually 
quoted as being relative to some cutoff momentum, often IGeV. 
However analysis of the GCM |10| gives rise to a different mass formula in which 

p = 2A J ^es{q')c{Q')crs{q') (32) 

with the factor c{s) = B{s)'^ /{sA[s)'^ + B(s)'^), and the expansion of the m dependence of 
B{q) from (12) 

B{q) Irn^o +m = B{q) \m=Q +es{q)m + 0{m^) (33) 

defines the current mass enhancement factor es{q). Despite the apparent difference be- 
tween (31) and (32) Langfeld and Kettner |Q have shown that they are indeed equal. 



9. Conclusions 

We have summarised here some of the facets of the GCM and its relationship to the fun- 
damental theory (QCD) and to various phenomenological models which have been prof- 
itably employed in modelling QCD. Significantly the GCM permits not only the derivation 
of these phenomenologies, but also particular expressions for the numerous parameters 
that arise in these models, such as the MIT bag constant, its modification within nuclear 
matter, and the ChPT parameters. Fundamentally the GCM is non-local and so the pre- 
dictions do not involve divergences and the consequent renormalisations, once the gluon 
propagator is given its scale dependence. 

The authors thank Professors K. Saito and A.W. Thomas and Dr A.G. Williams for the 
organisation of this workshop. 
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Figure 2 Variation of MIT Bag constant B{a)/B(l) wrt a. Partial restoration of chiral 
symmetry is described by cr < 1. 
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Table 1: Hadronic Observables 



Observable 


ml 

Theory 


Expt. /Theory 


Fitted observables 






u 


93.00MeV 


93.00MeV 


ai meson mass 


1230MeV 


1230MeV 


TT meson mass 


138.5MeV 


138.5MeV 


K meson mass (for mg only) 


496MeV 


496MeV 


Predicted observables 






{mu + md)/2 


6.5MeV 


6.0MeV 


rus 


135MeV 


130MeV 


uj meson mass 


804MeV 


782MeV 


Cq vr — vr scattering length 


0.1634 


0.21 ± 0.01 


Og TT — TT scattering length 


-0.0466 


-0.040 ± 0.003 


aj TT — vr scattering length 


0.0358 


0.038 ± 0.003 


0° TT — TT scattering length 


0.0017 


0.0017 ± 0.003 


vr — vr scattering length 


-0.0005 


not measured 


pion charge radius 


0.55fm 


0.66fm 


nucleon-core mass 


1390MeV 


~1300MeV 


constituent quark rms size 


0.59fm 


- 


chiral quark constituent mass 


270MeV 


- 


u/d quark constituent mass 


SOOMeV 


- 340MeV 


s quark constituent mass 


525MeV 


~510MeV 


0"*" diquark rms size 


O.TSfm 




0"*" diquark constituent mass 


692MeV 


>400MeV 


I"'' diquark constituent mass 


1022MeV 




0~ diquark constituent mass 


1079MeV 




1" diquark constituent mass 


1369MeV 




MIT core bag-constant 


(154MeV)4 


(146MeV)'^ 


MIT nucleon-core mass (no cm corr.) 


1500MeV 


~1300MeV 
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